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Purpose: To investigate the association with ocular biometric parameters in myopia-associated single nucleotide 
polymorphisms (SNPs) of the gap junction protein delta 2 (GJD2), insulin-like growth factor- 1 (IGF1) and hepatocyte 
growth factor (HGF) genes in two geographically different Chinese cohorts. 

Methods: In 814 unrelated Han Chinese individuals aged above 50 years including 362 inland residents and 432 island 
dwellers, comprehensive ophthalmic examinations were performed. Three SNPs, including GJD2 rs634990, IGF I rs62 14, 
and i/GFrs3735520, were genotyped. Genetic association with ocular biometric parameters was analyzed in individual 
cohorts, using linear regression controlled for sex and age. Common associations shared by the two cohorts were revealed 
by meta-analysis. 

Results: Meta-analysis showed that GJD2 rs634990 alone was not associated with any biometric parameters (adjusted 
p>0.645). The T allele of IGF 7 rs6214 was specifically associated with thicker lens (P±SE=0.055±0.022, adjusted 
p=0.034). The A allele of HGF rs3735520 was associated with longer vitreous chamber depth (P±SE=0. 143+0.060, 
adjusted p=0.050). Significant interaction between HGF rs3735520 and GJD2 rs634990 was found in association with 
axial length and vitreous chamber depth (adjusted p=0.003 and 0.033, respectively), and possibly with spherical error 
(adjusted p=0.056). 

Conclusions: Our endophenotyping analysis showed differential association between selected myopia-associated genes 
and ocular biometric parameters in our Chinese cohorts, which may underline substantial but diversified effects of these 
genes and their interaction on the development of eye structure and etiology of myopia. 



Myopia is one of the most common causes of visual 
impairment [ 1 -5] . It is estimated that about 33 . 1 % of the USA 
population is affected by this disorder [5]. The prevalence of 
myopia in China has been reported to be even higher, and up 
to 80% of Chinese children can have myopia [6,7]. Severe 
myopia is often linked to clinical complications [8], even 
permanent visual loss [9] . Myopia is an ultimate manifestation 
resulting from changes of eye structure or compartment in the 
optical path, which consists of cornea, anterior chamber, lens, 
and vitreous chamber [10,11]. These biometric parameters 
and the myopia disorder itself have been shown to have large 
genetic predisposition, implicating that these genetic 
determinants of ocular parameters can possibly influence the 
risk to myopia by controlling ocular development [12-15]. 
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Recently genome-wide association studies (GWAS) on 
quantitative traits have been successfully identified gene and 
variants associated with myopia. Variants at chromosome 
15ql4 and 15q25 have been reported to be associated with 
myopia and refractive error in two independent Caucasian 
GWAS [16,17]. Among these variants, single nucleotide 
polymorphisms (SNPs) of the gap junction protein delta 2 
(GJD2) gene at 15ql4 was reported to be more significantly 
associated with high myopia compared to SNPs at 15q25 in 
Japanese [18]. The GJD2 gene encodes connexin 36, a 36 kDa 
protein, which is a member of the connexin gene family and 
is highly expressed in mouse and human retina [19]. The 
connexin family can possibly be involved in ocular 
development and various eye diseases [20]. The quantitative 
trait association of GJD2 with refractive error thus remains to 
be investigated in Chinese. In addition to connexins, growth 
factors also play a substantial role in ocular development, and 
may influence biometric parameters [21]. The insulin-like 
growth factor 1 (IGF I) gene within the myopia 3 (high grade, 
autosomal dominant, MYP3) locus [22], has been reported to 
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be associated with myopia in Caucasians [23]. Expression of 
IGF1 mRNA in chicken ocular tissues can be affected by 
myopic or hyperopic defocus [24]. Likewise, association of 
the hepatocyte growth factor gene (HGF) with myopia has 
also been reported in Chinese [25] and Caucasians [26]. But 
quantitative trait association of both growth factor genes with 
ocular biometric parameters has not yet been studied. It 
remains to be investigated whether these three myopia- 
associated genes affect ocular development. 

In the current study we investigated the association of 
three myopia-associated genes, GJD2, IGF1, and HGF, and 
their interaction with eye biometric parameters in two Chinese 
cohorts. Our findings may suggest the substantial role of these 
genetic polymorphisms in shaping eye structure and 
development of myopia. 

METHODS 

Patient recruitment and clinical information: This study was 
approved by the Ethics Committee of Joint Shantou 
International Eye Center, Shantou, China and was conducted 
in accordance with the Declaration of Helsinki. Written 
consent was obtained from each participating subject after 
explanation of the nature of the study. 

The study subjects included 814 unrelated Han Chinese 
living all their lives in two geographical different regions in 
Southeastern China: 362 unrelated inland residents aged 50 
and older, recruited from senile cataract surgical patients at 
Joint Shantou International Eye Center in Shantou (STM), and 
432 unrelated local dwellers aged 50 and older, recruited from 
Nan'ao Island (NAI). The eyes with the following conditions 
were excluded: any history or symptom of Marfan' s 
syndrome, ocular trauma, ocular surgery, macular epiretinal 
membrane, macular edema, macular hemorrhage, glaucoma, 
or retinal detachment. Eye biometric parameters were 
documented for all study subjects. 

All participants received comprehensive ophthalmic 
examination including best-corrected visual acuity, slit-lamp 
biomicroscopy of anterior segment and retina with mydriasis, 
Refractive parameters including astigmatism, corneal 
curvature, spherical error, and cylindrical error were 
measured by auto refractometer (RK-F1 Refractometer/ 
Keratometer; Canon, Inc., Tochigi, Japan). Spherical 
equivalent was calculated as spherical error plus half of 
cylindrical error. Astigmatism was calculated as the 
difference between anterior and posterior cornea curvatures, 
and corneal curvature was calculated as the mean of the two. 
Axial length, anterior chamber depth, lens thickness, and 
vitreous chamber depth measured by A-scan ultrasound 
biometry (ODM 2200; Tianjin Maida Medical Technology 
Co., Ltd., Tianjin, China). The central corneal thickness 
(CCT) was measured ultrasonically (IOPac 20Mhz 
Pachymeter; Heidelberg Engineering, Heidelberg, Germany). 
Eyes with prior surgical history or low data quality were 
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excluded. For 557 individuals with bilateral data available, the 
means of biometric parameters was used to represent the data 
from both eyes. For 257 individuals with data from OD or OS 
data unavailable, data from the contralateral eye were used. 
Peripheral blood was collected from all participants, and 
genomic DNA was extracted using the Qiamp Blood Kit 
(Qiagen, Hilden, Germany). 

SNP genotyping: Three SNPs including rs634990 in GJD2, 
rs6214 in IGF-1, and rs3735520 in HGF were genotyped by 
Taqman SNP Genotyping assay (Applied Biosystems, Inc. 
[ABI], Foster City, CA) following the protocol suggested by 
the manufacturer. 

Statistical analysis: Hardy- Weinberg Test of each SNP was 
conducted using Haploview version 4.2 [27]. Gender 
difference between the two cohorts was compared using x 2 
tests, and age and biometric parameters were compared using 
non-parametric Mann-Whitney U test. Linear regression 
implemented by the R statistical language version 2.12.1 was 
used to analyze quantitative trait association for each 
individual cohort separately, controlling gender and age as 
described in previous studies [28,29]. The additive genetic 
model was used, assuming a trend per copy of the minor allele. 
Homozygous major, heterozygous, and homozygous minor 
genotypes were coded as 0, 1 , and 2 in the regression. Effect 
size±standard error (p±SE) of per copy of minor allele was 
calculated for each SNP accordingly. With the homozygous 
major genotypes as the reference (0 x P), heterozygous and 
homozygous minor genotypes were estimated to account for 
lxp and 2xp changes of biometric parameters, respectively. 
To identify common associations shared by the two Chinese 
cohorts, meta-analysis was further performed using fixed- 
effect models and inverse variance weighting methods 
implemented by METAL [30]. Bonferroni correction for 
multiple comparisons was applied to adjust meta-analysis p- 
values. 

RESULTS 

Demographic and clinical data: The distribution of refractive 
parameters and axial ocular dimensions in both STM and NAI 
cohorts were shown in Figure 1 and Figure 2. As Summarized 
in Table 1, comparison between the two Chinese cohorts 
showed significantly lower female proportion and older mean 
age in STM. It also revealed significant difference in both 
refractive parameters and axial ocular dimensions (all Mann- 
Whitney U test p<0.044). The STM cohort was in average 
more myopic with longer mean axial length, anterior chamber 
depth, and vitreous chamber depth, and thicker central cornea. 

Single gene association: None of the SNPs genotyped in the 
current study showed deviation from Hardy-Weinberg 
Equilibrium in either STM or NAI cohort (all p-value >0.05), 
and thus were subsequently included in further association 
study. The three SNPs showed similar minor allele 
frequencies between the two Chinese cohorts (42.6%— 49.0%, 
Table 2). 
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Figure 1. Distribution of refractive parameters in both the inland (STM) and island (NAI) cohorts. Histogram of the STM cohort is shown in 
light blue and that of the NAI cohort is in semitransparent red. 



Additive genetic models assuming a trend per copy of the 
minor allele were first used to test the association between 
biometric parameters and genotypes in each gene alone by 
using both eye data. As shown in Table 3 and Table 4, 
quantitative association analysis showed that GJD2 



rs634990 was association with central corneal thickness ((3 
±SE=-9.386±3.517 urn, p=0.008) in cohort STM. The 
association was not consistent in cohort NAI ((3 
±SE=0.819±2.131 jum, p=0.701), and became insignificant in 
meta-analysis (adjusted p=0.965). GJD2 rs634990 was not 
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Figure 2. Distribution of axial ocular dimensions in both both the inland (STM) and island (NAI) cohorts. Histogram of the STM cohort is 
shown in light blue and that of the NAI cohort is in semitransparent red. 



associated with any other refractive parameter or ocular 
dimension (all meta-analysis adjusted p>0.645). 

For IGF 1 rs6214, association between its minor allele T 
and corneal curvature was detected in STM (p±SE=0.23±0. 1 1 
D, p=0.029, Table 3). The association was insignificant in 
cohort NAI (p=0.629) and meta-analysis (adjusted p=0.1 1 1). 
The same allele T of IGF1 rs6214 showed a trend of 
association with longer lens thickness in both STM and NAI 



cohorts (p±SE=0.049±0.041 mm, p=0.240; and (3 
±SE=0.06±0.03 mm, p=0.027, respectively, Table 4). The 
association remained significant in meta-analysis ((3 
±SE=0.055±0.022 mm, adjusted p=0.034). No significant 
effects of IGF1 rs6214 was found on any other biometric 
parameters (all p>0.05). 

For 7/GFrs3735520, its minor allele A showed effects of 
negative spherical error and spherical equivalent, and longer 
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Table 1. Demographic information and clinical features of the study subjects. 



Category 


> 1 1V1 


IV AT 


P 


Gender 








Male 




1 1 A 

119 


^- A A A 1 

< 0.001 


Female 


ZZJ 


j j j 




Age (Year) 








Mean 


/l.O 


/CO ^ 

62.3 


^ A A A 1 

< 0.001 


(MJ) 




(y.2) 




Spherical error (D) 








Mean 


a a 

-0.9 


A "5 

-0.3 


A AA 1 

0.001 


(MJ) 


(3.5) 


(2.3) 




Cylindrical error (D) 








Mean 


A A 

—0.4 


A 1 

—(J. 3 


A A/I /I 

0.044 




(\ 1 ^ 
(1.1) 


(U.9J 




Spherical equivalent (D) 








Mean 


—1.2 


A /I 

-0.4 


^ r\ AA1 

< 0.001 




(4.1) 


(2.3) 




Astigmatism (D) 








Mean 


1 A 

— 1.0 


A *7 

—o. / 


^ r\ AA1 

< 0.001 




(o.y) 


(l.O) 




Curvature (D) 








Mean 


/I /I o 

44.3 


/I A A 

44.0 


A AA/T 

0.006 


(MJ) 


(1.!)) 


(!.:>) 




Axial length (mm) 








Mean 


23. o 


22. / 


< 0.001 


(MJ) 


(i./j 


(i.i) 




Central corneal thickness (jwm) 








Mean 


C A A H 

544.6 


531.2 


^ r\ r\(\ 1 

< 0.001 


(MJ) 




(30. /J 




Anterior chamber depth (mm) 








Mean 


3.2 


2.6 


< 0.001 


(SD) 


(0.4) 


(0.3) 




Lens thickness (mm) 








Mean 


4.4 


4.5 


0.004 


(SD) 


(0.6) 


(0.4) 




Vitreous chamber depth (mm) 








Mean 


16.2 


15.7 


< 0.001 


(SD) 


(1.6) 


(1.0) 





* x 2 tests were used for gender ratio comparison, and Mann-Whitney U tests were used for comparison of age and biometric 
parameter between the two cohorts. 



axial length in STM (p±SE=-1.03±0.40 D, p=0.011; (3 
±SE=-1.14±0.43 D, p=0.009 and (3±SE=0.35±0.13 mm, 
p=0.006, respectively, Table 3). These associations did not 
reach statistical significance in cohort NAI (p=0.298, 0.402 
and, 0.831, respectively). Meta-analysis did not found 
significance in these associations (adjusted p=0.147, 0.198, 
and 0.342, respectively). The same A allele of rs3735520 A 
showed a trend of association with longer vitreous chamber 
depth both in STM (p±SE=0.46±0. 14 mm, p=0.001), and NAI 
(p±SE=0.075±0.069 mm, p=0.277). And the association was 
marginally significant in meta-analysis ((3 



±SE=0.148±0.062 mm, adjusted p=0.050). No association of 
HGF rs3735520 was found with any other biometric 
parameters (all p>0.05). 

The same analysis was also performed using one eye data 
(Appendix 1 and Appendix 2), and the findings were 
comparable to the results above using both eye data. 

Gene-gene interaction: As shown in Table 5 and Table 6, 
meta-analysis of two-locus interaction was performed for the 
association of the three genes with ocular biometric 
parameters. By using meta-analysis, significant interaction 
between GJD2 rs634990 and HGF rs3735520 was revealed 
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in association with axial length and vitreous chamber depth 
(p±SE=-0.298±0.090, adjusted p=0.003 and (3 
±SE=-0.223±0.088, adjusted p=0.033, respectively). With the 
interaction item included in the full linear regression model, 
HGF rs3735520 showed significant effects on axial length 
and vitreous chamber depth (p±SE=0.373±0.104, adjusted 
p=0.001 and p±SE=0.359±0.103, adjusted p=0.001, 
respectively), and GJD2 rs634990 was associated with axial 
length (p±SE=0.231±0.096, adjusted p=0.049). Marginal 
significant interaction between GJD2 rs634990 and HGF 
rs3735520 was also found in association with spherical error 
and spherical equivalent (p±SE=0.540±0.231, adjusted 
p=0.056 and p±SE=0.559±0.250, adjusted p=0.075, 
respectively). When the interaction item included in the full 
linear regression model, HGF rs3735520 showed significant 
effects on axial length and vitreous chamber depth (p 
±SE=-0.804±0.226, adjusted p=0.006 and p 
±SE=-0.828±0.285, adjusted p=0.011, respectively). No 
significant interaction of IGF1 rs6214 with either of the other 
two genes was found in meta-analysis. 

DISCUSSION 

In the current study involving two geographically different 
Chinese cohorts, our results showed suggestive association of 
IGF I with lens thickness, and HGF with vitreous chamber 
depth. Hence our findings provided new insight into the roles 
of these myopia-associated genes in the development of 
different eye components in our Chinese cohorts and possibly 
in etiology of related eye diseases such as myopia. 

Ocular development and myopia can be shaped by 
genetic and environmental factors [31,32]. In the current 
study, dramatic differences in baselines of ocular biometric 
parameters were found between an inland cohort STM and an 
island cohort NAI. Such difference could be due to lower 
female proportion and older mean age in STM. Moreover, it 
could be due to variable environmental exposure and lifestyle 
between the two. In contrast, the two Chinese cohorts have 
close genetic background in the three genes investigated in the 
current study. In spite of such difference in trait baselines 
between the two cohorts, common genetic correlation with 
ocular biometric parameters could be detected within each 
individual cohort, and further confirmed by a meta-analysis 
approach. These findings suggested that intrinsic genetic 
factors contributed to variations of ocular biometric 
parameters that could not be explained by environmental 
factors. 

The quantitative trait association studies have been used 
to delineate genetic predisposition in these disease-related 
biometric parameters. Previously Solouki et al. [17] reported 
chromosome 15ql4 spanning SNP rs634990 in GJD2 to show 
genome- wide significance for association with refractive 
error in a Dutch population-based GWAS. The C allele of 
rs634990 was recently reported to confer risk to myopia in 
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Japanese [18]. Although our cohorts showed similar minor 
allele frequency of rs634990 compared to the Hapmap Han 
Chinese data and the Japanese cohort, its association with 
spherical equivalent or other refractive parameters was not 
detected. The current findings might indicate ethnic difference 
in genetic predisposition of myopia between our Chinese 
cohorts and other reported populations. Moreover, our two- 
locus analysis results implicated that GJD2 could play a role 
in myopia etiology by interacting with other myopia- 
associated genes in ocular development and association with 
biometric parameters. 

The genotype frequencies of rs6214 in our Chinese 
cohorts were similar to the reported Han Chinese of 
Hapmap data. IGF1 rs6214 was specifically associated with 
lens thickness in our Chinese cohorts. The minor allele of 
IGF1 rs6214 was correlated with 0.07 mm increase of lens 
thickness in our meta-analysis Chinese cohort, which account 
for about 1.56 D change in refractive error according to 
previously reported approximately 0.045 mm/D change in 
lens thickness [33,34]. The lens of adult human accounts for 
about one third of the total refractive power in the eye [35]. 
Although correlation of IGF1 with refractive error was not 
detected, the change of lens thickness could still potentially 
affect the ultimate refractive error. Previously, rs6214 was 
reported to be associated with both high myopia and myopia 
in an international Caucasian cohort [23]. Animal studies have 
implicated the role of the IGF1 in lens development. IGF1 has 
previously been reported to induce lens cell elongation and 
specialized crystallin gene expression in embryonic chicken 
eyes [36]. The association of IGF 7 with lens thickness but not 
with other ocular dimensions, constellated with the existing 
genetic association of IGF1 with myopia, possibly implicated 
its specific role in refractive myopia. 

In contrast to IGF 7, HGF 'was specifically associated with 
axial length and vitreous chamber, but not lens thickness in 
our meta-analysis. The minor allele A was correlated with 
increased chamber depth in the meta-analysis. Axial length is 
one of the major determinants of refractive error, and accounts 
about 50% variance of spherical equivalent [15]. Vitreous 
chamber is the largest compartment in the eye, and its depth 
accounts for the largest proportion of axial length. These 
findings could explain the previous report of HGF as a high 
myopia-associated gene in the Chinese population [25]. 
Intriguingly, HGF exhibited significant interaction with 
another myopia-associated gene GJD2, which also 
contributed to the genetic association with axial length and 
vitreous chamber depth. Notably such interactive effects were 
also implicated in spherical error and spherical equivalent but 
not cylindrical error, and HGF was significantly associated 
with these two parameters when interaction was considered. 
HGF probably interact with GJD2 to control the axial 
dimension and thus influence refractive parameters, which 
possibly explain its association with myopia. Axial length 
change has been estimated to be 0.35 mm/D in myopia [37], 
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and thus the effect size of 0.377 mm per copy of rs3735520 
minor allele was expected to account for approximately 1.07 
D change of refractive error, which was close to the observed 
value of 0.828 D. And homozygous minor genotypes of 
rs with 2 copies of minor alleles could account for 

1.656 D change of refractive error. Interestingly, in the 
Nepalese population HGF was recently reported to be 
associated with primary angle-closure glaucoma [38], in 
which patients were usually featured by shorter axial length 
and vitreous chamber depth. The SNP rs3735520 was 
associated with serum HGF level in normal individuals [39], 
suggesting its possible function link to gene expression. Taken 
together, HGF is probably involved in development of the 
posterior eye segment, and consequently in spherical error and 
axial myopia. 

Myopia is characterized by major clinical features 
including negative refractive error and elongated eye axial 
length. However, both of these two features are ultimate 
phenotypes depending on various genes modulating the 
anatomic development of the eye. The differential correlation 
of myopia-associated genes with refractive error and axial 
ocular dimensions in the current study thus underlined the 
importance of endophenotyping in myopia genetics study. 
Firstly different genes or gene sets could be responsible for 
specific endophenotypes. Moreover, genes that controlled 
axial length could be of special interest. It has been reported 
that these genes account for approximately 50% of the 
variation in spherical equivalence [15]. Secondly, our data 
further pointed to a substantial role of interaction between 
these genes such as HGF and GJD, in genetic studies of 
myopia endophenotypes. 

In the current study, we reported differential phenotype- 
genotype correlations between myopia-associated genes and 
eye biometric parameters in the Chinese population. IGF1 was 
associated with lens thickness, HGF was associated with 
vitreous chamber depth, and the interaction between HGF and 
GJD2 was associated with axial length, vitreous chamber 
depth and possibly spherical error. These findings provided 
new information in the diversified functional role of these 
susceptibility genes in myopia etiology and ocular 
development. 

ACKNOWLEDGMENTS 

We would like to acknowledge that this study was in part 
supported by the research grants from Natural Science 
Foundation of Guangdong Province, China (No. 
8 1 5 1 503 1 020000 1 9), National Natural Science Foundation of 
China (No. 81000397), Science and Technology Planning 
Project of Guangdong Province, China (No. 
20 10B03 1600130) and Joint Shantou International Eye 
Center, Shantou University/The Chinese University of Hong 
Kong (No. 08-001, 08-002, 09-014, 10-020, 10-021, and 10- 
022). We would also like to thank Ms. Pancy Tarn, Dr. Lijia 
Chen and Dr. Gary Yam for helpful discussions. 



© 2012 Molecular Vision 

REFERENCES 

1 . Bourne RR, Dineen BP, Ali SM, Noorul Huq DM, Johnson GJ. 

Prevalence of refractive error in Bangladeshi adults: results 
of the National Blindness and Low Vision Survey of 
Bangladesh. Ophthalmology 2004; 111:1150-60. [PMID: 
15177965] 

2. Dandona R, Dandona L, Srinivas M, Giridhar P, McCarty CA, 

Rao GN. Population-based assessment of refractive error in 
India: the Andhra Pradesh eye disease study. Clin Experiment 
Ophthalmol 2002; 30:84-93. [PMID: 11886410] 

3. Kempen JH, Mitchell P, Lee KE, Tielsch JM, Broman AT, 

Taylor HR, Ikram MK, Congdon NG, O'Colmain BJ. The 
prevalence of refractive errors among adults in the United 
States, Western Europe, and Australia. Arch Ophthalmol 
2004; 122:495-505. [PMID: 15078666] 

4. Sawada A, Tomidokoro A, Araie M, Iwase A, Yamamoto T. 

Refractive errors in an elderly Japanese population: the Tajimi 
study. Ophthalmology 2008; 115:363-70e3. [PMID: 
18243904] 

5. Vitale S, Ellwein L, Cotch MF, Ferris FL 3rd, Sperduto R. 

Prevalence of refractive error in the United States, 1999— 
2004. Arch Ophthalmol 2008; 126:1111-9. [PMID: 
18695106] 

6. He M, Zheng Y, Xiang F. Prevalence of myopia in urban and 

rural children in mainland China. Optom Vis Sci 2009; 
86:40-4. [PMID: 19104465] 

7. Tang WC, Yap MK, Yip SP. A review of current approaches to 

identifying human genes involved in myopia. Clin Exp 
Optom 2008; 91:4-22. [PMID: 18045248] 

8. McBrien NA, Gentle A. Role of the sclera in the development 

and pathological complications of myopia. Prog Retin Eye 
Res 2003; 22:307-38. [PMID: 12852489] 

9. Saw SM. How blinding is pathological myopia? Br J 

Ophthalmol 2006; 90:525-6. [PMID: 16622078] 

10. Curtin BJ, Karlin DB. Axial length measurements and fundus 

changes of the myopic eye. I. The posterior fundus. Trans Am 
Ophthalmol Soc 1970; 68:312-34. [PMID: 5524211] 

11. Curtin BJ, Karlin DB. Axial length measurements and fundus 

changes of the myopic eye. Am J Ophthalmol 1971; 
71:42-53. [PMID: 5099937] 

12. Ashton GC. Segregation analysis of ocular refraction and 

myopia. HumHered 1985; 35:232-9. [PMID: 4029963] 

13. Lyhne N, Sjolie AK, Kyvik KO, Green A. The importance of 

genes and environment for ocular refraction and its 
determiners: a population based study among 20^15 year old 
twins. Br J Ophthalmol 2001; 85:1470-6. [PMID: 11734523] 

14. Vitart V, Bencic G, Hayward C, Herman JS, Huffman J, 

Campbell S, Bucan K, Zgaga L, Kolcic I, Polasek O, 
Campbell H, Wright A, Vatavuk Z, Rudan I. Heritabilities of 
ocular biometrical traits in two Croatian isolates with extended 
pedigrees. Invest Ophthalmol Vis Sci 2010; 51:737-43. 
[PMID: 19875653] 

15. Dirani M, Shekar SN, Baird PN. Evidence of shared genes in 

refraction and axial length: the Genes in Myopia (GEM) twin 
study. Invest Ophthalmol Vis Sci 2008; 49:4336-9. [PMID: 
18552384] 

16. Hysi PG, Young TL, Mackey DA, Andrew T, Fernandez- 

Medarde A, Solouki AM, Hewitt AW, Macgregor S, 
Vingerling JR, Li YJ, Ikram MK, Fai LY, Sham PC, Manyes 



776 



Molecular Vision 2012; 18:765-778 <http://www.molvis.org/molvis/vl8/a82> 

L, Porteros A, Lopes MC, Carbonaro F, Fahy SJ, Martin NG, 
van Duijn CM, Spector TD, Rahi JS, Santos E, Klaver CC, 
Hammond CJ. A genome-wide association study for myopia 
and refractive error identifies a susceptibility locus at 15q25. 
Nat Genet 2010; 42:902-5. [PMID: 20835236] 

17. Solouki AM, Verhoeven VJ, van Duijn CM, Verkerk AJ, Ikram 

MK, Hysi PG, Despriet DD, van Koolwijk LM, Ho L, Ramdas 
WD, Czudowska M, Kuijpers RW, Amin N, Struchalin M, 
Aulchenko YS, van Rij G, Riemslag FC, Young TL, Mackey 
DA, Spector TD, Gorgels TG, Willemse-Assink JJ, Isaacs A, 
Kramer R, Swagemakers SM, Bergen AA, van Oosterhout 
AA, Oostra BA, Rivadeneira F, Uitterlinden AG, Hofman A, 
de Jong PT, Hammond CJ, Vingerling JR, Klaver CC. A 
genome-wide association study identifies a susceptibility 
locus for refractive errors and myopia at 15ql4. Nat Genet 
2010; 42:897-901. [PMID: 20835239] 

18. Hayashi H, Yamashiro K, Nakanishi H, Nakata I, Kurashige Y, 

Tsujikawa A, Moriyama M, Ohno-Matsui K, Mochizuki M, 
Ozaki M, Yamada R, Matsuda F, Yoshimura N. Association 
of 15ql4 and 15q25 with high myopia in Japanese. Invest 
Ophthalmol Vis Sci 2011; 52:4853-8. [PMID: 21436269] 

19. Sohl G, Joussen A, Kociok N, Willecke K. Expression of 

connexin genes in the human retina. BMC Ophthalmol 2010; 
10:27. [PMID: 20979653] 

20. White TW, Bruzzone R. Intercellular communication in the eye: 

clarifying the need for connexin diversity. Brain Res Brain 
Res Rev 2000; 32:130-7. [PMID: 10751662] 

21. Klenkler B, Sheardown H. Growth factors in the anterior 

segment: role in tissue maintenance, wound healing and 
ocular pathology. Exp Eye Res 2004; 79:677-88. [PMID: 
15500826] 

22. Young TL, Ronan SM, Alvear AB, Wildenberg SC, Oetting 

WS, Atwood LD, Wilkin DJ, King RA. A second locus for 
familial high myopia maps to chromosome 12q. Am J Hum 
Genet 1998; 63:1419-24. [PMID: 9792869] 

23. Metlapally R, Ki CS, Li YJ, Tran-Viet KN, Abbott D, Malecaze 

F, Calvas P, Mackey DA, Rosenberg T, Paget S, Guggenheim 
JA, Young TL. Genetic association of insulin-like growth 
factor- 1 polymorphisms with high-grade myopia in an 
international family cohort. Invest Ophthalmol Vis Sci 2010; 
51:4476-9. [PMID: 20435602] 

24. Penha AM, Schaeffel F, Feldkaemper M. Insulin, insulin-like 

growth factor-1, insulin receptor, and insulin-like growth 
factor- 1 receptor expression in the chick eye and their 
regulation with imposed myopic or hyperopic defocus. Mol 
Vis 2011; 17:1436-48. [PMID: 21655358] 

25. Han W, Yap MK, Wang J, Yip SP. Family-based association 

analysis of hepatocyte growth factor (HGF) gene 
polymorphisms in high myopia. Invest Ophthalmol Vis Sci 
2006; 47:2291-9. [PMID: 16723436] 

26. Yanovitch T, Li YJ, Metlapally R, Abbott D, Viet KN, Young 

TL. Hepatocyte growth factor and myopia: genetic 



© 2012 Molecular Vision 

association analyses in a Caucasian population. Mol Vis 
2009; 15:1028-35. [PMID: 19471602] 

27. Barrett JC, Fry B, Mailer J, Daly MJ. Haploview: analysis and 

visualization of LD and haplotype maps. Bioinformatics 
2005; 21:263-5. [PMID: 15297300] 

28. Chen J, Tsang SY, Zhao CY, Pun FW, Yu Z, Mei L, Lo WS, 

Fang S, Liu H, Stober G, Xue H. GABRB2 in schizophrenia 
and bipolar disorder: disease association, gene expression and 
clinical correlations. Biochem Soc Trans 2009; 37:1415-8. 
[PMID: 19909288] 

29. Chen JH, Wang D, Huang C, Zheng Y, Chen H, Pang CP, Zhang 

M. Interactive effects of ATOH7 and RFTN1 in association 
with adult-onset primary open angle glaucoma. Invest 
Ophthalmol Vis Sci. 2012 [PMID: 22222511] 

30. Wilier CJ, Li Y, Abecasis GR. METAL: fast and efficient meta- 

analysis of genomewide association scans. Bioinformatics 
2010; 26:2190-1. [PMID: 20616382] 

3 1 . Wilson A, Woo G. A review of the prevalence and causes of 

myopia. Singapore Med J 1989; 30:479-84. [PMID: 2694377] 

32. Hornbeak DM, Young TL. Myopia genetics: a review of current 

research and emerging trends. Curr Opin Ophthalmol 2009; 
20:356-62. [PMID: 19587595] 

33. Koretz JF, Cook CA, Kaufman PL. Accommodation and 

presbyopia in the human eye. Changes in the anterior segment 
and crystalline lens with focus. Invest Ophthalmol Vis Sci 
1997; 38:569-78. [PMID: 9071209] 

34. Dubbelman M, Van der Heijde GL, Weeber HA. Change in 

shape of the aging human crystalline lens with 
accommodation. Vision Res 2005; 45:117-32. [PMID: 
15571742] 

35. Gordon RA, Donzis PB. Refractive development of the human 

eye. Arch Ophthalmol 1985; 103:785-9. [PMID: 4004614] 

36. Bassnett S, Beebe DC. Localization of insulin-like growth 

factor-1 binding sites in the embryonic chicken eye. Invest 
Ophthalmol Vis Sci 1990; 31:1637-43. [PMID: 2167301] 

37. Atchison DA, Jones CE, Schmid KL, Pritchard N, Pope JM, 

Strugnell WE, Riley RA. Eye shape in emmetropia and 
myopia. Invest Ophthalmol Vis Sci 2004; 45:3380-6. [PMID: 
15452039] 

38. Awadalla MS, Thapa SS, Burdon KP, Hewitt AW, Craig JE. 

The association of hepatocyte growth factor (HGF) gene with 
primary angle closure glaucoma in the Nepalese population. 
Mol Vis 2011; 17:2248-54. [PMID: 21897747] 

39. Burdon KP, Macgregor S, Bykhovskaya Y, Javadiyan S, Li X, 

Laurie KJ, Muszynska D, Lindsay R, Lechner J, Haritunians 
T, Henders AK, Dash D, Siscovick D, Anand S, Aldave A, 
Coster DJ, Szczotka-Flynn L, Mills RA, Iyengar SK, Taylor 
KD, Phillips T, Montgomery GW, Rotter JI, Hewitt AW, 
Sharma S, Rabinowitz YS, Willoughby C, Craig JE. 
Association of Polymorphisms in the Hepatocyte Growth 
Factor Gene Promoter with Keratoconus. Invest Ophthalmol 
Vis Sci 2011; 52:8514-9. [PMID: 22003120] 



777 



Molecular Vision 2012; 18:765-778 <http://www.molvis.org/molvis/vl8/a82> 



©2012 Molecular Vision 



Appendix 1. 



Association of GJD2 rs634990, IGF I rs6214, and HGF 
rs3735520 with refractive parameters using one eye data. To 
access the data, click or select the words "Appendix 1." This 



will initiate the download of a compressed (pdf) archive that 
contains the file. 



Appendix 2. 



Association of GJD2 rs634990, IGF I rs6214, and HGF 
rs3735520 with axial ocular dimensions using one eye data. 
To access the data, click or select the words "Appendix 2." 



This will initiate the download of a compressed (pdf) archive 
that contains the file. 
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